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Abstract: The NEXT experiment aims to observe the neutrinoless double beta decay 
of 186 Xe in a high-pressure xenon gas TPC using electroluminescence (EL) to amplify the 
signal from ionization. One of the main advantages of this technology is the possibility to 
reconstruct the topology of events with energies close to Qpp. This paper presents the first 
demonstration that the topology provides extra handles to reject background events using 
data obtained with the NEXT-DEMO prototype. 

Single electrons resulting from the interactions of 22 Na 1275 keV gammas and electron- 
positron pairs produced by conversions of gammas from the 228 Th decay chain were used 
to represent the background and the signal in a double beta decay. These data were used 
to develop algorithms for the reconstruction of tracks and the identification of the energy 
deposited at the end-points, providing an extra background rejection factor of 24.3 ±1.4 
(stat.)%, while maintaining an efficiency of 66.7 ± 1. % for signal events. 
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1 Introduction 


The observation of neutrino oscillations by several experiments proves that neutrinos have 
mass, as described, for instance, in Ref. [1]. This, however, motivates a number of questions 
about the nature of neutrinos and their relation to the other fermions: How can this mass be 
accommodated in the Standard Model? What is the absolute scale of neutrino mass? Does 
CP violation take place in neutrino oscillations? Could the neutrino be its own antiparticle? 

All other fermions in the Standard Model (quarks and charged leptons) have their mass 
described by a Dirac mass term alone. Neutrino mass could also be generated by the same 
mechanism, however, the neutrino, being chargeless, can also have a Majorana mass term. 
Majorana neutrinos would be indistinguishable from their antiparticle and may lead to 
processes violating lepton number. Such processes may be responsible for the generation 
of the observed matter-antimatter asymmetry of the Universe via leptogenesis [2], The 
Majorana mass of the neutrino also provides an elegant explanation of the smallness of 
neutrino mass via the seesaw mechanism [3]. 

Currently, the most sensitive experimental method to establish that neutrinos are 
Majorana particles is the search for neutrinoless double beta decay (/3/30^). This is a 
hypothetical, very rare nuclear transition in which a nucleus with Z protons decays into a 
nucleus with Z+2 protons and the same mass number, A, emitting two electrons that carry 
essentially all the energy released (Qpp)- While the two-neutrino mode of the double beta 
decay has already been measured in a number of isotopes, the zero-neutrino mode remains 
unobserved. 

The experimental signature of a neutrinoless double beta decay are two electrons with 
total kinetic energy equal to Qpp- For this reason, to be able to make the measurement, an 
experiment must be optimised simultaneously for energy resolution and the rejection of 
non-signal events with energies similar to Qpp- Energy resolution is key not only to limit the 
number of events caused by natural radiation which enter the region of interest but to limit 
the pollution from the intrinsic background caused by the two neutrino mode. The signature 
of a (3/32is event only differs from signal in the sum of the energies of the two electrons 
and, therefore, improved separation is only made possible by improved energy resolution. 
Backgrounds from ambient sources can be further reduced through a combination of fiducial 
cuts and particle identification. 

The NEXT experiment seeks to make a first measurement of /3/30U in 136 Xe using a high 
pressure gas Time Projection Chamber (TPC) with electroluminescent (EL) read-out. It 
has been designed to provide good energy and spatial resolution to identify separated tracks 
and increased ionization (‘blobs’) at their ends. To this end, it uses two different planes for 
energy measurement and tracking. Photomultipliers behind the cathode detect the primary 
scintillation light and allow the energy measurement by detecting the electroluminescent 
light. An array of silicon photomultipliers (SiPMs) behind the anode provides the spatial 
information for the topological analysis of the events using the EL light. The f3f30v search 
will be carried out using NEXT-100 which will contain 100 kg of 136 Xe gas at 15 bar and 
is described in detail in Refs. [4, 5]. The first phase of the experiment, called NEW, and 
deploying 10 kg (and 20% of the sensors), is currently being commissioned at the Laboratorio 
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Figure 1. A f3/30v event (left) and a single electron background event from a 2.44 MeV 214 Bi 
gamma (right) in Monte Carlo simulation (both events simulated at 15 bar gas pressure). 

Subterraneo de Canfranc. The principle of operation of NEXT-100 and NEW and the 
necessary know-how has been developed using NEXT-DEMO, a large scale prototype which 
operated at the Instituto de Fisica Corpuscular in Valencia (Spain) with ~1.5 kg of natural 
xenon at a pressure of 10 bar (for a detailed description of the prototype, see Refs. [6, 7]). 

The use of a topological particle identification based on the expected signature of a 
double electron (signal) event compared to that of a single electron (background) produced 
by the interaction of high energy gammas is presented here. Using the Monte Carlo 
simulation of the NEXT-DEMO prototype and data taken with NEXT-DEMO a first 
demonstration of the power of the method has been made. This involved the comparison of 
single electron tracks originating from the photoelectric interaction of 22 Na gammas and 
double electron tracks from the pair production of the 2.614 MeV gamma from 208 T1. 

The paper is organized as follows. The topological signal and reconstruction algorithms 
are described in Sec. 2. The data analysis is described in Sec. 3 and the results obtained 
are presented and discussed in Sec. 4. The paper ends with the conclusion in Sec. 5. 

2 Topological signature in NEXT 

Electrons (and positrons) moving through xenon gas lose energy at an approximately fixed 
rate until they become non-relativistic. At the end of the trajectory the 1/v 2 rise of the 
energy loss (where v is the speed of the particle) leads to a significant energy deposition 
in a compact region, which will be referred to as a ‘blob’. The two electrons produced 
in double beta decay events appear as a single continuous trajectory with a blob at each 
end (Fig. 1-left). Background events from single electrons, however, typically leave a single 
continuous track with only one blob (Fig. 1 -right). The use of this topological signature to 
eliminate background in f3/3 experiments was pioneered by the Caltech-Neuchatel-PSI 
Collaboration in the Gotthard Underground Laboratory [8], using a gaseous 136 Xe TPC 
with multiwire read-out, with a fiducial mass of 3.3 kg of 136 Xe at a pressure of 5 atm. 
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Figure 2. Distribution of the energy deposition in the two blob candidates of a track, for Monte 
Carlo /3/30v (left) and 208 T1 (right) events, simulated with the NEXT final detector [9]. Blob 
candidates are defined as the sum of all charge within a 2 cm radius of each end-point of the tracks. 

At the 15 bar pressure of NEXT, the two electrons emitted in a neutrinoless double beta 
decay tend to leave a single track of about 15 cm length. Energy deposition is approximately 
constant along the length of the track other than at the two extremes where the electrons 
deposit more than 20% of the event energy, split between both blobs. The main background 
in NEXT comes from high energy gammas emitted in 208 T1 and 211 Bi decays, which occur 
naturally in the detector materials as part of the 232 Th and 238 U chains, entering the active 
volume of the detector. These gammas convert in the gas through photoelectric, Compton 
and pair production processes producing, typically, more than one electron separated in 
space, the tracks of which have just one blob of energy at one extreme. Figure 2 shows 
the distribution of the energy deposited in the blob candidates of the tracks for signal 
and background events, simulated in NEXT-100. As can be seen, the energies of the blob 
candidates have closer values for signal events with the blob candidate with less energy at 
much lower energy for background. 

In a previous publication [7], only the average position of the track was reconstructed, 
using the barycentre calculated with the total integrated signal in each SiPM. The work 
presented here uses a more sophisticated algorithm to reconstruct the topology left in the 
detector by a charged particle. 

2.1 Reconstruction of tracks in an EL TPC 

Reconstruction of tracks in an electroluminescent TPC is complicated not only by the 
diffusion of the charge cloud during drift but also by the nature of the read-out. Scintillation 
light is produced over the whole width of the EL gap (5 mm in NEXT-DEMO) spreading 
the signal from a single electron over a time inversely proportional to the drift velocity 
within the gap (~ 2 /rs). Additionally, the EL light is produced isotropically and, therefore, 
the signal produced by the passage of an electron through the gap is expected to arrive at 
the tracking plane (~7.5 mm behind the anode) over the area defined by the intersection of 
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the plane with the sphere of light. 

In a previous paper [10], the NEXT collaboration demonstrated that a ‘point-like’ 
deposition of charge due to the absorption of a K a X-ray is expected to be detected over a 
transverse region which can be parameterised as a two dimensional Gaussian with a standard 
deviation of ~8 mm where the spread due to EL light production is the dominant effect 
with subdominant contributions from transverse diffusion of the charge. Longitudinally, 
the expected spread has a noticeable dependence on the drift distance since the diffusion 
dominates. K a events are expected to have widths in z with standard deviations of between 
0.5 mm, for very short drifts, and 1.7 mm at the drift field settings used here. In order to 
optimise the reconstruction of tracks these values must be taken into account by dividing the 
signal information into appropriate time slices and using charge information from clustered 
SiPM channels. 

2.2 Hit-finder algorithm 

Once a signal-like pulse is selected using the energy plane, the corresponding information for 
tracking must be analyzed in the SiPM plane. The SiPM noise is mainly due to electronic 
noise in the laboratory, which is higher than the dark current. The distribution of noise per 
time sample has been measured for each SiPM through dedicated runs without external 
signal and used to eliminate from the analysis those time slices which exhibit charge below 
a certain noise level. 

The total charge in a time integrated section (slice) of the tracking plane is required to 
be above a minimum threshold in order to be considered by the algorithm. This requirement 
has the purpose of eliminating the residual noise, due to fluctuations, without affecting the 
signal of low energy depositions in the detector. Typically, slices with total charge below 
threshold are located at the beginning and end of the pulse. In these slices, the position of 
the sensors with charge is unrelated to the energy depositions reconstructed in the other 
slices of the pulse. In Fig. 3 -left, an example of this kind of slice is shown. After this first 
requirement, the information for each time slice is passed to the hit-finder algorithm. 

The algorithm used in this work searches for clusters around local maxima and then 
proceeds iteratively, selecting first the channel with maximum charge and forming a cluster 
with the first ring of sensors around it. If, and only if, the ring is fully occupied (taking 
into account any dead channels or edge effects) and each individual SiPM exhibits a charge 
higher than a threshold, the cluster information is used to build a hit, whose x and y 
position are reconstructed as the barycentre of the charge information. The information 
already used is then removed and the procedure repeated until no SiPM above threshold are 
left. The need for this threshold can be understood by considering a typical slice located in 
the centre of a pulse, like the one shown in Fig. 3 -right. The light seen at the level of a 
few photoelectrons across the whole SiPM array has undergone multiple reflections before 
arriving at the sensors and, therefore, provides no relevant position information. The energy 
recorded for the corresponding time slice in the energy plane is then divided between all 
hits found according to the proportion of charge in each hit. This charge is then corrected 
according to the detector calibration described in Ref. [10]. 
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Figure 3. Tracking plane signals for two arbitrary pulse slices, as an illustration of the hit-finder 
algorithm. (Left) A slice that does not pass the minimum energy cut, (right) a slice with clear signal 
(see text for more details). 


2.3 Voxelization and track identification 

Once a set of hits is found, a connectivity criterium must be defined so that the hits 
belonging to each separate particle can be grouped into tracks. The procedure is as follows: 
first, the active volume is divided into 3D pixels, known as ‘voxels’, of fixed dimensions. 
Each voxel is given an energy equal to the sum of the energies of all the hits which fall 
within its boundaries. The collection of voxels obtained in such a way can be regarded 
as a graph, defined as a set of nodes and links that connect pairs of nodes. Two voxels 
can then be considered connected if they share a face, an edge or a corner, with each pair 
of connected voxels being given a weight equal to the geometric distance between their 
centres. Next, the “Breadth First Search” (BFS) algorithm (see, for instance, Ref. [11]) 
is used to group the voxels into tracks and to find their end-points and length. The BFS 
algorithm is a graph search algorithm which finds the minimum path between two connected 
nodes, starting from one node and exploring all its neighbours first, then the second level 
neighbours and so on, until it reaches the second node. The BFS algorithm divides the 
voxels into connected sets, known as tracks and finds their end-points, defined as the pair of 
voxels with largest distance between them, where the distance of two voxels is the shortest 
path that connects them. The distance between the end-points is the length of the track. 

The choice of the voxel size is a compromise between a fine granularity and conservation 
of connectivity, which depends on the hit-finder algorithm in use. In this work, the best 
performance has been found for voxels of 1 x 1 x 1 cm 3 . Improvements in the hit-finder 
algorithm will allow for smaller voxels or different connectivity criteria and a more accurate 
reconstruction of tracks will be possible. However, one has to take into account that the 
current transverse and longitudinal diffusion in NEXT-DEMO, as well as transverse spread 
due to the EL production, affect the position of the ionization electrons during the drift, 
blurring the track, thus spatial resolution better than the typical diffusion does not provide 
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sufficient additional information for the additional complexity and cost to be worthwhile. 


3 Analysis 

The analysis was performed using data from two different calibration sources to study the 
topology of background-like and signal-like events in NEXT-DEMO. The active volume of 
the detector has a drift length of 30 cm and a hexagonal cross section of 8-cm apothem. 
Background-like events were studied using a l-/iCi 22 Na calibration source. The de-excitation 
of the first excited state of the 22 Na daughter isotope, 22 Ne, produces a gamma with an 
energy of 1.275 MeV. These gammas can produce a photo-electron which leaves a track in 
NEXT-DEMO of approximately 7 cm length (10 bar pressure). This track is contained in 
the TPC and is long enough to perform topological studies. 

Signal-like events were obtained using a 228 Th source. The 228 Th decay chain includes 
208 T1, the daughter of which, 208 Pb, is created in an excited state which de-excites emitting 
a 2.614 MeV gamma. This gamma can produce an electron-positron pair with a signature 
which mimics, except for the total energy, the topology of a f3(30v event. The positron 
produces a blob equivalent in energy to that of the electron, and, when it annihilates, emits 
two back-to-back 511 keV gammas. Due to the size and the pressure of NEXT-DEMO, 
there is a high probability for both gammas to escape the active volume. In this case, the 
energy deposited in the chamber is 1.592 MeV and the track left by the electron-positron 
pair is around 6 cm long. Although the energy of these signal-like events is higher than that 
of background-like events, the track is slightly shorter on average because both the electron 
and the positron, individually, have a lower energy than the photo-electron of the 22 Na 
decay. These two sources produce tracks of comparable size and energy hence providing 
two data sets with good characteristics for a topological study. 

3.1 Data set 

The 22 Na and 208 T1 data were taken under the following detector conditions: 10 bar pressure, 
667 V cm -1 drift field and 2.4 kV cm -1 bar -1 EL field. The sources were placed outside 
the lateral port positioned midway between cathode and anode. 

The 22 Na data were taken in November 2013 using a trigger on the scintillation signal 
optimised for the higher energy gamma compared to that described in [6]. An integrated 
exposure of ~84.4 hours was used in the analysis presented. The 208 T1 data were collected 
in May 2014 using a trigger on the EL signal. An event was required to contain at least 
one ionization signal which registered at least 21 000 photo-electrons (pe) in the central 
PMT. An integrated exposure of ~165.68 hours was used in this case. 

Dedicated Monte Carlo (MC) simulation samples were generated using NEXUS, a 
Geant4 [12, 13] based package developed by the NEXT collaboration that simulates the 
entire signal generation process. Samples of interacting gammas at the energy of that from 
22 Na and that from 208 T1 have been used for event selection studies. In order to study a 
possible source of background coming from cosmic muons, an additional sample of muon 
events with an energy of 4 GeV and an angular distribution following ~ cos 2 ($), according 
to Ref. [14], where 6 is the zenith angle, has been analyzed. 
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3.2 Event selection 

Charge pulses are identified as originating from primary scintillation (SI) or EL (S2) based 
on the criteria that they have a duration in time of between 1 and 3 ys for SI (>3 /is for 
S2), and that they have integrated charge of ^0.5 pe (^10 pe for S2) in the average of all 
the PMTs. In this way both the higher energy events of interest for the analysis and the 
xenon X-rays are included. X-ray pulses allow for the energy scale to be easily validated 
between separate runs. Only events with one, and only one, S 1-like signal are accepted as 
potential signal to avoid event pile up and for ease of calculation of the 2 coordinate. Before 
the data are passed to the hit-finder algorithm, the pulses identified in the events must first 
pass a number of pre-selection requirements designed to remove any spurious charge due 
to instrumental effects. These include requirements on the drift to ensure that no S2 has 
deposited charge within 2 cm of either the cathode or anode; on the relative mis between 
the recorded signal in the PMTs rejecting S2s with PMT rms greater than 20% of the 
mean signal; on the ratio between anode and cathode charge; and a loose transverse fiducial 
requirement on the integrated SiPM charge (|x| and \y\ ^ 60 mm). Table 1 shows the 
effect of this pre-selection on the data and Monte Carlo samples. The table also shows the 
reductions in the total data samples for the two sets. The pre-selection requirements remove 
far more data events since they are designed to remove events that are not induced by the 
source which are not present in the MC sample, such as interactions in the cathode buffer 
region and small sparks. However, it can be seen that at the level of the final requirement 
the samples are at a comparable level for the 22 Na data with the difference seen for the 
228 Th sample explained by the presence of single electron events in the region of interest 
(ROI) as described in Sec. 4.2, validating the filtering method. 

All events which passed the preselection were then presented to the hit-finder algorithm 
described in Sec. 2.2 and subjected to a final fiducial requirement based on the position of 
the reconstructed hits. No event may have any hits reconstructed outside the cylinder with 
cross section defined by R = \J x 2 + y 2 = 50 mm and z limits as defined by the fiducial cut 
described above (20 < z < 280 mm). This requirement ensures that entering backgrounds 
and events only partially contained are removed from the data set. All remaining simulated 
cosmic muon events were rejected by this fiducial requirement. 




Remaining sample proportion 


Selection 


MC 

Data 


208 ^ 

to 

to 

p 

208 qq 

22 Na 

Pre-selection 

0.451T0.002 

0.0769T0.0001 

0.329T0.001 

0.0608T0.0002 

xy inside fid. veto 

0.244T0.002 

0.411T0.001 

0.286T0.001 

0.316T0.002 

Energy in the ROI 

0.470T0.005 

0.0071T0.0002 

0.274T0.003 

0.063T0.002 

One track identified 

0.955T0.003 

0.976T0.005 

0.887± 0.004 

0.847T0.010 

Two blobs found 

0.658T0.007 

0.219T0.013 

0.558± 0.006 

0.243T0.014 


Table 1 . Event selection efficiencies for each cut requirement described in the text evaluated for 
signal MC plus muon sample and for data. Errors are statistical. 
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After the final fiducial requirement the energy of each S2 was reconstructed using the 
position dependent calibration described in Ref. [10]. Events were then voxelised and tracks 
were constructed using the algorithms described in Sec. 2.3. Using voxels of 1 x 1 x 1 cm 3 , 
track sections were considered in three dimensions and, when two reconstructed tracks of 
at least four voxels were separated by at most one voxel, they were merged to form a single 
track. Figure 4 shows the energy spectra for both sources for all events with one merged 
track (details on the one-track requirement will be discussed later in the text). In these 
figures, the photoelectric peak of the interaction of the 22 Na 1.275 MeV gamma and the 
double escape peak resulting from pair production induced by the 2.614 MeV gamma can 
be seen. The 208 T1 data were taken 3 months after the calibration run resulting in a poorer 
energy resolution. In NEW and NEXT-100 regular calibration runs will be taken to avoid 
similar decalibration. 

In order to exclude events which are not photoelectric interactions (in 22 Na data) or 
pair production events (in the case of 208 Tl) a ROI was selected for each dataset. The 
22 Na ROI was chosen by fitting a Gaussian to the 1275 keV peak and considering all 
events with energy within ±3 a of the mean obtained from the fit. A similar selection 
was made for 208 T1 data taking into account that in the region of the double escape peak 
single electron events from Compton scattering are also expected. The fit was performed 
parameterising the fall of the Compton edge with an exponential and the signal with a 
Gaussian distribution. The fit predicts that 25.7 ± 0.6% of events are single electrons and 
74.3 ± 0.6% are pair production events in a region within ±3 a of the mean obtained from 
the fit. The exponential background model fits well in the regions immediately before and 
after the peak, and the variation of the range used in the combined fit predicts errors on 
the signal to background ratio of less than 1%. 




Energy (keV) 


Figure 4. Energy spectra for the 22 Na source (left) and the 228 T1 i source (right), for one-track only 
events passing the selection criteria described in the text. 


Those events which have tracks with energies within the ROI are then filtered according 
to the number of tracks identified. Events with more than one track are not selected for 
the final sample. At the energies considered here both signal and background are expected 
to produce only one track due to which this requirement has little effect. However, in the 
f5/30u analysis described in Ref. [9] selecting only one track proves important and is kept 
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here for completeness. As can be seen in Table 1, 89% of the 208 T1 double escape peak 
candidates and 85% of the 22 Na 1275 keV gamma candidates fulfill this requirement. 

For each track in the ROI, the blobs were sought at their end-points. A blob candidate 
is defined as a group of voxels which fall within a set radius from a track end-point. Figure 
5 shows, for both datasets, the energies of the blob candidate with lower energy plotted 
against that of the other blob candidate. It can be seen that the single-electron-dominated 
sample from 22 Na tends to have less energy in its lower energy blob candidate. A slight 
shift in the energies of the higher energy blob candidates is seen also, which is expected to 
be due to calibration effects. 


>1200 

o> 



200 400 600 800 1000 1200 

Higher energy blob cand. energy (keV) 



200 400 600 800 1000 1200 

Higher energy blob cand. energy (keV) 


Figure 5. Energy distribution at the end-points of the tracks coming from 22 Na decay (left) and 
those coming from the 208 T1 decay (right) for 2 cm radius blob candidates. 


4 Results 

Single electrons are only expected to have a high energy deposition at one end of their track 
as opposed to double electrons which are expected to have large energy depositions at each 
end of the track. An example of single electron track candidate is shown in Fig. 6 and a 
typical example of double electron track candidate is shown in Fig. 7. Figure 5 shows that 
a requirement on the energy of the less energetic blob candidate is the simplest criterium 
to separate signal ( 208 T1 e + e _ ) from background ( 22 Na e _ ) events. Optimisation of this 
requirement uses the standard figure-of-merit e/Vb, where e is the signal efficiency and b is 
the fraction of background events that survives the requirement in the signal region. 

4.1 Monte Carlo studies 

The evaluation of the figure of merit was performed as a function of the blob candidate 
radius and the minimum energy required for the low energy blob candidate, using MC 
data (see Fig. 8 -left). An additional consideration is the probability of blob candidate 
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Figure 6. Example of a 22 Na gamma track. The zx and zy projections (on top) and the yx 
projection in the lower frame are drawn. The energy scale is in keV (color coded) and the distances 
are in mm. The single electron candidate exhibits one blob at one end of the track. 


region overlap. In Fig. 8- right, the fraction of events with overlapping blob candidates is 
shown for different radii, for both 22 Na and 208 T1 samples. For larger radii, blob candidate 
overlap occurs in a higher percentage of events, particularly in the case of signal events 
which tend to be shorter. A radius of 2 cm is found to be optimal if one considers that 
it is large enough to contain the energy deposition at the end of the tracks while being 
small enough, compared to the average track length, to limit the number of tracks in which 
voxels are considered in both blob candidates to ~3%. For blob candidates of 2 cm radius, 
an optimum requirement on the minimum blob candidate energy is found at 210 keV, as 
can be seen in Fig. 8 -left. From the range of values that give very similar figures of merit, 
that with higher signal efficiency was selected. 

Applying this requirement to the minimum blob candidate energy allows 65.9% of the 
remaining 208 T1 events to pass and 21.9% of the 22 Na events. Figure 9 shows the efficiency 
for acceptance/survival of the two blob requirement for the 208 T1 and 22 Na samples in 
the ROI as a function of the observed energy. Events in the 22 Na sample are expected to 
consist of single electrons and the survival efficiency is not expected to be a strong function 
of the electron energy in the small range of the ROI, see Fig. 9 -left. The final sample 
of 208 T1 events is expected to be a sum of the pair production and single electrons from 
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Figure 7. Example of a double escape peak candidate from the 208 T1 data sample. The zx and zy 
projections (on top) and the yx projection in the lower frame are drawn. The energy scale is in keV 
(color coded) and the distances are in mm. The double electron candidate exhibits two blobs, one at 
each end of the track. 




Figure 8. Left: the figure of merit, e/y/b as a function of the blob candidate radius and the 
minimum blob candidate energy. Right: proportion of events with blob candidate overlapping, as 
a function of the blob candidate radius. The step-like behaviour that appears in both plots is an 
artificial consequence of the 1-cm size of voxels. For certain values of radius, the number of voxels 
inside the blob candidate, thus the proportion of events with overlap, does not change. 
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Compton scattering with the composition dependent on the event energy, therefore it is 
expected that the acceptance efficiency should exhibit energy dependence, as illustrated in 
Fig. 9- right. The efficiency of the two blob requirement on these two event categories was 
studied using separated Monte Carlo event samples, as shown in Fig. 10. As expected, the 
relative population of these categories varies strongly with the event energy, but the effect 
of the requirement for a given event category does not vary significantly within the ROI. 
The survival factor for the simulated single electron events in the 208 T1 sample is found to 
be consistent with that measured for the 22 Na sample, between 20% — 30%, whereas the 
efficiency of acceptance of double track events is at the level of 65% — 70%. 



Energy (keV) 



Figure 9. Two-blob cut efficiency for the Monte Carlo 22 Na (left) and 208 T1 (right) samples. The 
histograms represent the energy of the event, while the points are the efficiency of the two-blob cut 
for the events in the bin. 



Energy (keV) 



Figure 10. Two-blob cut efficiency for the Monte Carlo 208 T1 sample, selecting events with pair 
production only (left) and events without pair production (right). The histograms represent the 
energy of the event, while the points are the efficiency of the two-blob cut for the events in the bin. 


Validation of the MC and reconstruction methods has been performed. The MC was 
found to reproduce the topological features found in data to a high degree of accuracy 
as can be seen in Fig. 11 -left, for the track lengths and in Fig. 11-right for the energy of 
the higher energy blob candidate shown for data and MC. In Fig. 11 -right we see a slight 
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systematic shift between MC and data, which we will study in greater detail in the future. 




Figure 11. Comparison between data and Monte Carlo simulations, for the reconstructed length of 
the tracks (left) and the reconstructed energy inside the blob candidate with higher energy (right). 
22 Na samples are used. Bins 12 and 13 in the bottom left plot are not shown due to lack of MC 
statistics (larger samples will be used for the future analysis). 


4.2 NEXT-DEMO results 

Applying the same selection to the data taken with NEXT-DEMO results in a proportion 
of 55.8 ± 0.9 (stat.)% of the 228 Th dataset passing the two-blob requirement and 24.3 ± 1.4 
(stat.)% for the 22 Na sample. In the case of 208 T1 this result does not represent an efficiency 
of identification of double escape peak events since the ROI contains both single electron 
tracks (Compton events) and double electron tracks (pair production events). Assuming 
the proportions of each interaction type predicted by the parameterisation of the spectrum 
described in Sec. 3.2 and the single-electron survival of the 22 Na data, an estimation of the 
signal efficiency can be made according to the following equality: 

^total ^bkg * /bkg /, .. 

^signal = -7-“-“ , l 4 - 1 ) 

/signal 

where e s i gna i and £bk g are the proportions of signal and single-electron events that pass 
the cut, etotai is the total proportion of events in the final sample, and /bkg and / s i gna i 
are the fractions of background and signal-like events in the sample subjected to the cut 
(/bkg +/signal = 1)- This procedure predicts that 66.7±0.9 (stat.) ± 0.3 (fit)% of the double 
escape peak events pass the final cut. The same procedure applied to the MC samples gives 
a result of 68.6 ± 0.8 (stat.)% efficiency for the double escape peak events. 

Fig. 12 -left shows the final two-blob selection for both data and MC, displaying signal 
efficiency for the double escape peak events and background rejection (defined as the 
proportion of events that does not pass the two-blob selection) varying the minimum energy 
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required for the lower energy blob. Fig. 12-right shows signal efficiency as a function of the 
value of the energy requirement. Data and MC simulations are in agreement to within 2 
statistical sigmas for all points, a further validation of the NEXT Monte Carlo. 




Figure 12. (Left) Signal efficiency as a function of the required minimum energy of the lower energy 
blob candidate. Both data and Monte Carlo simulation are shown and the values corresponding to 
the cut used in this analysis are highlighted. (Right) Signal efficiency as a function of background 
rejection (proportion of background events removed from the sample by the two-blob cut) varying 
the required minimum energy of the lower energy blob candidate. 


4.3 Discussion 

The power of the topological signature in rejection of background and acceptance of signal 
events in a high-pressure xenon gas TPC (NEXT-DEMO) is well described by the detailed 
simulation and reconstruction of the Monte Carlo events (the agreement between data and 
MC is better than 3%). The NEXT-DEMO MC samples have also been simulated using the 
‘fast simulation’ which forms voxels using a Gaussian smearing of the energies of the true 
deposits. The efficiencies obtained for the two-blob cut are in agreement within 5% with 
the results of this paper. This agreement allows for an informed comparison between these 
results and those of Ref. [9] where the fast simulation of NEXT-100 was used to analyse 
events from PPOv, 214 Bi, and 208 T1. 

Compared to the present analysis, the fast simulation and analysis of NEXT-100 
predict significantly improved background survival rate of 10% as opposed to 24% for the 
NEXT-DEMO configuration, at the same signal efficiency. A major difference which affects 
these results is the relative size of the two detectors. NEXT-100 has a drift length of 1.3 m 
with a circular cross section of ~1 m and will operate at 15 bar pressure making it easily 
large enough to contain electron tracks at energies similar to Qpp regardless of topology 
or orientation. NEXT-DEMO, on the other hand, is significantly smaller resulting in the 
fiducial cuts favouring more tortuous tracks which do not displace as far from their origin 
and are more difficult to reconstruct and more prone to blob candidate overlap. This bias, 
which will not be present in NEXT-100, is expected to account for the differences observed. 

The topological signatures and background rejection capabilities of the algorithm 
presented in this paper demonstrate the adequate performance for the proposed NEXT-100 
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experiment. Future experiments may require larger, one tonne scale, detectors and also 
better background rejection. The high potential of the topological signatures has been 
demonstrated by the Gotthard experiment, where only 3.5% of the background events were 
allowed to pass the final selection requirements, albeit with the detector configuration which 
compromised the energy resolution [8]. It is expected that future improvements in the 
design of the high-pressure xenon gas TPC and more accurate reconstruction algorithms 
will provide significant improvements above the level expected for NEXT-100 detector. 

5 Conclusions 

The possibility to identify single and double electron topologies in a xenon gas TPC has 
been demonstrated using the NEXT-DEMO detector. This topological separation provides 
a powerful handle for the rejection of the mainly single electron backgrounds faced by 
gaseous xenon (3(30v experiments. In the NEXT-100 experiment a reduction of an order of 
magnitude in background is expected using this method. 

Modelling signal using topologies left by electron-positron pairs induced by gammas 
from 208 T1 decay and background using single electrons induced by the de-excitation gamma 
produced as part of the 22 Na chain, the topological separation has been optimised and 
applied to both Monte Carlo and data. Considering only the final requirement where a 
minimum end-point energy of 210 keV was required, a signal efficiency of 66.7 ± 0.9 (stat.) 
±0.3 (fit)% was measured. The same selection results in 24.3 ± 1.4 (stat.)% of background 
events entering the final sample. A parallel study using the simulation of f3f30u events 
and backgrounds at that energy in NEXT-100 has found a similar efficiency, for roughly 
half of the background [9]. This difference can be accounted for considering the difference 
in size both of the events themselves and of the detectors. Higher energy events tend to 
be longer, making it easier to differentiate MIP-like sections from the end-points of the 
events improving the rejection of backgrounds at Qpp. The small radius of NEXT-DEMO 
coupled with the requirement of full containment also favours less extended events where 
the reconstruction is more complicated. The analysis will be repeated using data from the 
next stage of the NEXT experiment which uses radiopure materials in its construction and 
has a larger fiducial volume. 
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